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We investigate within the Gross-Pitaevski (GP) theory the formation of 

fringe patterns between two Bose condensates stored in a double well potential 

modelling the forces applied to the experimental systems. In the case of 

repulsive interactions between the atoms, we report the onset of interference 

structures similar to those observed in the experiments, after the release of the 

confining potential. Conversely, attractive interactions lead to the collapse of 

the condensate when the number of particles is larger than a critical value. We 

show that a bias field introduced in the center of the trap allows the system 

to avoid the blow-up of the density and gives rise to a periodic behavior of 

growth and decay of spatial modulations reminiscent of the Fermi-Pasta-Ulam 

recurrence. 



The successful experimental realizations ]T]-|J of Bose-Einstein condensation (BEC) in 
dilute alkali atomic vapors have sparkled a renewed interest in the study of weakly coupled 
Bose gases. At temperatures close to zero, these magnetically trapped Bose condensed clouds 
are well described by the condensate wavefunction ^(r, £) which obeys the GP equation 
initially introduced for superfluid helium [|5] 

^ = [ -^ y2 + c/| * |2 + ^ (r)] *' (1) 

with the normalization condition 



J \^(r,t)\ 2 dr = N, (2) 

where TV is the number of particles in the condensate. The interactions between atoms of 
mass m is described by the potential U. Since the thermal wavelength is larger than the 
characteristic length associated with binary collisions, this potential can simply be charac- 
terized by the s- wave scattering length: U = 4:7rh 2 a/m. The interaction is repulsive (a > 0) 
for 87 Rb and 23 Na and attractive (a < 0) for 7 Li. The confining potential is given by V ext (r). 
Many works have been devoted to the characterization of the ground state properties and 
collective excitations of such condensed gases confined by a single harmonic potential f^-Q . 

Recently interference fringes have been observed to form when two condensates of repul- 
sive atoms are allowed to expand and overlap UlQ| . Interference between bose condensates 



with a > has also been the subject of various theoretical studies fTT Hl3|1 - In this letter, we 
investigate the formation of such fringe patterns stressing the difference between the cases of 
attractive and repulsive interaction and we report new dynamical behaviors for bose gases 
with a < 0. In this framework, we numerically solve the GP equation to follow the time 
evolution of two condensates with an equal number of atoms initially stored in a double well 
potential. It consists of the superposition of a harmonic potential describing the magnetic 
trap and a gaussian barrier in the center modelling the laser sheet that cuts the condensate in 
half. As in the experiments, the iso-potentials of the harmonic trap are ellipsoidal with their 
long axis defining vertical direction z. The ratio Q = ^ < 1 of the corresponding angular 
frequencies of the trap in the vertical and transverse directions characterizes the anisotropy 
of the potential. The height of the barrier is chosen to make tunneling ineffective. This 
initial condition can thus be written in the form: 

^0, x, y, 0) = (j) L {z + d, x, y, 0) + (f> R (z - d, x, y, 0)e iff , (3) 

and S is the relative phase which proceeds from the breaking of the groupe symmetry and 
varies between different experimental runs. The 0's are ground state wavefunctions in each 
of the two asymmetric traps created by the laser. 
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We first consider the case of repulsive interaction. When TV is sufficiently large (Thomas- 
Fermi regime), the wavefunctions are then essentially determined by the balance between 
the particle interaction energy and the external potential as the contribution of the kinetic 
energy is negligible. The spatial extension L of (j) is much smaller than the initial sepa- 
ration 2d between the two condensates (Fig. la). When the double trap is switched off, 
the condensates expand ballistically and fringe patterns appear in the overlap region. This 
expansion is essentially dominated by the dispersive effects generated by the kinetic energy 
term that has been frustrated in the ground state (j) of the trap. Similar numerical expansion 
dynamics are indeed obtained for a = and a ^ 0. The wavelength of these interference 
structures increases linearly with time \(t) — ^ and the corresponding expansion velocity 
hn/md decreases for increasing values of the width of the gaussian barrier. The experimen- 
tal fringe period also becomes smaller for larger powers of the argon laser sheet and thus for 



larger initial distances between the condensates []10J . The essential dynamics can already be 
captured by ID numerical simulations along the z direction that produce the time develop- 
ment illustrated in Fig. la which presents strong analogies with the time-resolved diffraction 



patterns obtained in double-slit atomic experiments WM. These ID structures are described 



by the interference term I{z ) t) in the expression of the atomic density |^(r, t)\ 2 which takes 
the form: 

/occos(^ + 0). (4) 



Variation of the relative phase 5 only leads to a trivial shift of the whole structure ]TT[ . Taking 
advantage of the rotational symmetry about the vertical axis, two-dimensional simulations 
of Eq. (1), that corresponds for instance to a section along the O xz plane of the 3D system 



produce fringe patterns (Fig. lb) similar to those seen in the experiments on 23 Na [|T0| . 
These structures unambiguously demonstrate the coherence of those Bose condensate gas. 
The system chooses a phase for the macroscopic wavefunction in a process which is analogous 
to the appearance of coherent oscillations in a laser or a chemical reactor. Indeed the GP 
equation also corresponds to the strong dispersion limit of the complex Ginzburg-Landau 



equation which describes these symmetry breaking bifurcations [|T5|. The atomic cloud of 
5 10 5 atoms then behaves like a single entity to which the de Broglie relation can be applied. 

The dynamics is completely different in the case of attractive potentials U = — \U\ < 0. 
This arises because the uniform condensate wavefunction = 

exact solu- 
tion of the (GP) equation in the absence of a confining potential, may undergo a Benjamin- 
Feir instability ]16[ (n is the number density in the condensate). Owing to the form of the 



dispersion relation 



h imn \U\k 2 
Uh= 2^t R^"' (5) 



the state is unstable to long wavelength density perturbations with a wavenumber lying 
in the range 



2Jmn \U\ 

< k < V ' ' (6) 

It had therefore been claimed that BEC of attractive bosons was impossible as this 
condensation would have corresponded to a mechanically unstable state for which the com- 
pressibility would become negative fl?p . Nevertheless convincing experimental evidence for 



BEC in a gas of 7 Li atoms in a trap, consisting of permanent magnets ]TB[, has been reported 
recently @. In such a confined system, BEC can be achieved when the number of particles is 
sufficiently small to inhibit the destabilizing long wavelength perturbations. More precisely 
this occur when TV < N c = a(l/\a\) where I is the typical extension of the ground state in 
the trap and a is a number of order unity that is determined by the detailed form of the 
confining potential 0. When N > N c , the attractive potential overwhelms the zero point 
energy leading to the collapse of the wavefunction by which, in 2D and 3D, a singularity 
tends to form in a finite time [1^]2(J . A sufficient criterion for the collapse ( VPT criterion) 



also yields the critical number N c . In the collapsing cloud, new effects not contained in Eq. 
1 should be taken into account such as various intrinsic inelastic processes which lead to the 
heating of the sample. 

In order to study the interaction between two Bose condensates of attractive atoms, 



we consider the same configuration as preceedingly. Here however we maintain some con- 
finement at all times to preserve the integrity of the condensed phase. Initially two stable 
(N < N c ) independent condensates, located near the minima at z = ±d are separated by the 
Gaussian barrier. When this barrier is lowered a non-trivial dynamical behavior develops 
(Fig. 2). During the expansion along the z direction the unstable modes that were excluded 
from the small initial condensates are awakened by the increase of the size of the system. As 
a result, a structure appears that invades the whole confining region. Contrary to the case of 
repulsive interactions, the wavelength does not vary with time anymore and corresponds to 
that of the fastest growing mode in the unstable range (Eq. (6)). After reaching a maximum 
value, the amplitude of the modulated density decreases and the system returns to a state 
of localized condensates near the minima of the potential (Fig. 2). This process repeats 
periodically and its "superperiod" decreases when the intensity of the barrier is increased. 
This behavior is reminiscent of the well studied Fermi-Pasta-Ulam (FPU) recurrence and 



of other recurrences obtained, for instance, with the nonlinear Schrodinger equation pi 
that corresponds to the V ext = of Eq. (1). It is also analogous to the phenomenon of 
periodic alternation observed in nonlinear optics experiments |22| . Both phenomena consist 



of a periodically ordered sequence of quasi-stationary modes. A modal decomposition of 
the wave function on the most unstable modes reproduces the spatial periodic solutions of 
the GP equation. The superperiod which can be calculated analytically and expressed in 
terms of elliptic functions |23| is, as the wavelength, in good agreement with the numerical 



simulations. 

By moving the atoms away from the center along the z axis, the small potential barrier 
acts to lower the central density hump below the critical value N c and so prevents the ignition 
of the collapse of the global system since the cloud can be populated by more atoms before 
reaching the blow-up. The creation of a vortex in a rotating cloud can lead to a similar 
effect §. 

At larger values of N (but still N < A c ), more unstable modes are allowed, leading to 
the numerical observation of imperfect recurrences and of more complex behaviors involving 



a restricted range of wavenumbers [24j]. When eventually exceeds N c each condensate 
collapses separately. 

The main effect of the harmonic potential is to confine the system and consequently to 
reduce the number of unstable modes taking part in the recurrence. When the anisotropy 
of the confining potential is important, ft < 1, the fringes are aligned along the tranverse 
direction, and the structure is periodic in the vertical direction. When the system is confined 
by a weaker potential, transverse modes also fall into the unstable range of wavenumbers (Eq. 
(5)), leading to cellular structures with spots covering the whole confinement region. They 
compete with fringes during the recurrent motion, each one taking its turn at dominating 
the solution profile. 

Atomic clouds of near zero temperature, much smaller than the critical condensation 
temperature, may now be obtained experimentally. In this regime, dissipative effects and 
thermal fluctuations do not affect the macroscopic dynamics anymore. Furthermore, owing 
to the presence of the trapping potential, zero temperature fluctuations may be strongly 
suppressed |P5| . Therefore BEC of attractive bosons could provide an ideal candidate to 
experimentally observe FPU type recurrence phenomena. 

In the presence of an initial asymmetric repartition (n i and n^) of the atoms between the 
two condensates in the double trap potential, we further report numerical observations of an 
oscillatory exchange of atoms, analogous to the Josephson effect. Numerical simulations in 
both ID and 2D exhibit the already described FPU recurrence (Fig. 3), now accompanied 
by a periodic exchange of particles locked at half the recurrence frequency. This locking 
phenomena is observed independently of the initial density repartition. 

Contrary to the case of repulsive interactions for which the fringe patterns can be inter- 
preted as resulting from interferences between de Broglie waves, condensates of attractive 
atoms can exhibit for modulated initial conditions a recurrent behavior that is intimately 
related to the Benjamin- Feir instability of the uniform macroscopic wavefunction. 
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Figure 1. Interference patterns for two freely expanding Bose condensates of 23 N a atoms 
(a = 4.9 nm, N = 2 10 5 ) initially separated by d = 40 \im in a double well, 
(la) space-time plot of the density |\^| 2 obtained by numerically integrating the Gross- 
Pit aevski equation (Eq. 1) switching off the double- well at t — 0. Space is along the z 
direction and time runs downward up to t = 40ms. The linear increase of the wavelength 
with time is visible. 

(lb) density plot in the xz plane after 40 ms when the wavelength equals of 15 fim. The 
gray scale interpolates between maximum (dark) and minimum (light) densities. 

Figure 2. Recurrent behavior of a Bose condensate of 7 Li atoms (a = —1.45 nm ) N = 700) 
evolving in a double well potential: sequence of density plots (plane xz) spanning the re- 
currence " superperiod" . At time t — 0, two identical condensates were captured in their 
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respective well and separated by an adequate barrier. This barrier is lowered to produce the 
exhibited time evolution. 

Figure 3. Josephson-type behaviour for a BE condensate of 7 Li atoms (a = — 1.45nra, 
N = 700). Space-time plot of the density along the z axis during a single exchange of 
particles along the z axis when the initial transients have died out. 
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